arXiv:1506.05577vl [astro-ph.SR] 18Jun2015 


Peak Bagging of red giant stars observed by Kepler, first 
results with a new method based on Bayesian nested sampling 


Enrico Corsaro* ® and Joris De Ridder' 

Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium 


Abstract. The peak bagging analysis, namely the fitting and identification of single os¬ 
cillation modes in stars’ power spectra, coupled to the very high-quality light curves of 
red giant stars observed by Kepler, can play a crucial role for studying stellar oscillations 
of different flavor with an unprecedented level of detail. A thorough study of stellar oscil¬ 
lations would thus allow for deeper testing of stellar structure models and new insights in 
stellar evolution theory. However, peak bagging inferences are in general very challenging 
problems due to the large number of observed oscillation modes, hence of free parame¬ 
ters that can be involved in the fitting models. Efficiency and robustness in performing 
the analysis is what may be needed to proceed further. For this purpose, we developed a 
new code implementing the Nested Sampling Monte Carlo (NSMC) algorithm, a power¬ 
ful statistical method well suited for Bayesian analyses of complex problems. In this talk 
we show the peak bagging of a sample of high signal-to-noise red giant stars by exploit¬ 
ing recent Kepler datasets and a new criterion for the detection of an oscillation mode 
based on the computation of the Bayesian evidence. Preliminary results for frequencies 
and lifetimes for single oscillation modes, together with acoustic glitches, are therefore 
presented. 


1 Introduction 

Since the first detection of non-radial oscillations in cool giant stars Q, the asteroseismology of field 
red giants (RGs) has encountered a substantial increase, especially thanks to the advent of the space- 
based photometric missions CoRoT (e.g. rrsiTTi ) and Kepler (e.g. IfTTO ). where the latter allowed 
also for the study of RGs in open clusters (e.g. HSU). 

The discovery of so-called mixed modes d has led to a significant improvement in our under¬ 
standing of the internal structure and evolution of RGs. In particular, the characteristic spacing in 
period of mixed modes frequencies provides a direct way to disentangle H-shell and He-core burning 
RGs II2I20II and has been used already to classify about 13,000 targets observed in the Kepler field of 
view Ga. Therefore, RGs are among the most interesting and useful types of stars to refine models 
and test thoroughly codes of stellar structure and evolution. 

With the present 4-years long datasets available from Kepler for many RGs, we have the possibility 
to perform very detailed analyses aimed at measuring asteroseismic properties such as frequency, 
amplitude, and lifetime, of each individual oscillation mode observed in the Fourier spectra of the stars. 
This particular analysis, usually known as peak bagging plays a key role to exploit the full potential 
of the informative power spectra of the stars. However, due to the large number of modes populating 
the power spectrum of a RG (up to -100), the peak bagging usually turns into a computationally 
demanding analysis II9|15|10|4| . The complexity of the analysis is even higher considering that the 
mode identification of all the oscillation peaks is necessary for a proper modeling of the oscillations. 
This implies that also testing the significance of individual oscillation peaks becomes an important task 
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Fig. 1. Resulting background fit of the star KIC 12008916, as derived by DIAMONDS. The original PSD is shown 
in gray, while a smoothed version with boxcar width set to zly/5, with Av taken from m , is shown as a black line 
to guide the eye. The red thick line represents the background model without the Gaussian envelope, while the 
red dotted line accounts for the additional Gaussian component. The individual components of the background 
model are shown by blue dot-dashed lines. 


within the process. For this purpose, (H developed a new code termed DIAMONDS that implements the 
nested sampling Monte Carlo algorithm ED, a sampling algorithm well suited for model comparison 
and efficient high-dimensional inferences in a robust Bayesian perspective fH. 


2 Observations and Data 

To select the final sample of stars for the analysis we started from the original sample of RGs observed 
in long cadence (LC) by Kepler and studied by IfTTl and by O, consisting in a total of 1111 stars. 

From the original sample we only considered those stars with power spectra having an oscillation 
envelope with frequency of maximum power Vmax ^ 110/iHz, according to the stellar population 
study done by (161, with y^ax values provided by ifTlI . We then checked for available measurements 
of their g-mode period spacing AUi, as measured by (201. This cross-match led to 21 stars observed 
by Kepler from QO till Q17.1, resulting in a frequency resolution of Vbm - 0.008yuHz for each power 
spectral density (PSD). The Kepler light curves were corrected according to O. 

In this work we focus on a low-mass low-luminosity RGB star, KIC 12008916, a particular type 
of RG that is part of a small population of stars (about 5 % of the entire population of RGs observed 
by Kepler), having high signal-to-noise ratio (SNR) oscillations. Due to the clear oscillation pattern 
and the large number of modes, this star is well suited for testing stellar evolution theory and pulsation 
codes, as well as to constrain stellar formation rates in our galaxy. 


3 Peak Bagging model and results 

The first step for the peak bagging analysis of the stars is to estimate the background signal in their 
PSD. We follow the probabilistic study done on a large sample of stars by (121, who found that two 
separate granulation components with a characteristic frequency close to the that of the region contain¬ 
ing the oscillations, together with a third component related to long-term variations, best represent the 
background signal of the RGs selected. Thus, the background model that we adopt can be expressed as 
the sum of three super-Forentzian profiles, a flat noise component, and a Gaussian envelope to model 
the power excess of the oscillations. The resulting background fit for KIC 12008916 is shown in Fig.[^ 
The subsequent step in analyzing the stellar PSDs is to adopt an adequate fitting model for the 
oscillation pattern contained in the region of the power excess. In this work we adopt a method similar 
to that used by n (see Sect. 6.3). We restrict our frequency range of analysis to the region contain¬ 
ing the oscillations, which we identify as y^ax ± 3.5(Tenv, with (Xenv the standard deviation from the 
Gaussian envelope used to fit the background. We thus take into account some substantial differences 
characterizing the oscillation pattern in RGs, as explained in the following. 
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In the oscillation pattern of all the RGs, one can observe forests of modes between consecutive 
quadrupole and radial modes, which are known as dipole mixed modes, namely oscillations with a 
mixed character between that of pure p and g modes. Since g modes have the highest mode inertia 
because they propagate in high-density regions, their lifetime r is significantly longer than that of pure 
pressure modes (e.g. see 0). As a consequence, even with an observing time Tobs > 4 years now 
made available from Kepler, mixed modes with a g-dominated character appear still unresolved, while 
others having a more p mode-like character are partially or sometimes even fully resolved. 

In this work we distinguish the resolved or partially resolved peaks, for which Tobs ^ t, from those 
that are unresolved, for which Tobs ^ t. In the former case, the oscillation peak profile that we adopt 
is that of a Lorentzian and it is given by 


^res,0 (V) 


AlKnrp) 


( 1 ) 


where Aq, Tq, vq are the amplitude, the linewidth and the frequency, respectively, and represent the 
three free parameters to be estimated during the fitting process. According to the Fourier analysis, an 
oscillation peak that is not resolved has a profile represented as lO : 
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where Hq and vq are the height in PSD units and the central frequency in pRz of the oscillation peak, 
respectively, and must be estimated during the fitting process, while Vbm is the frequency resolution 
introduced in Sect. and it is a fixed value. In this case fitting the height is preferred since it is an 
observable and we have no information about the linewidth of the peak, while fitting a Lorentzian 
profile would represent an inadequate model as it involves a spurious fitting parameter represented by 
the linewidth. Then, the final peak bagging model can be represented as a mixture of Lorentzian and 
sinc^ functions. 

Some of the preliminary results for the peak bagging analysis of KIC 12008916 are shown in 
Fig.[^ where one can observe a depression of the linewidths of dipole mixed modes in correspondence 
to the position of the frequency of maximum power y^ax, and in Fig.[^ with the evidence of acoustic 
glitches for the quadrupole and radial modes. 


3.1 Peak significance 

Testing the significance of an oscillation peak in the PSD of a star is a crucial aspect that has to be 
considered for providing a reliable set of modes that can be used for the modeling of the oscillations, 
hence for investigating the stellar structure and the evolution of the star. For this purpose, following 
i) (Sect. 6.5) we consider a detection probability for a single peak, p^, defined as 


Pb = 


Sa + Sb 


( 3 ) 


where 6 a and 6 b are the Bayesian evidences for the models excluding and including, respectively, 
the oscillation peak to be tested, and are computed with Diamonds. The detection limit is usually set 
to Pb > 0.99 as derived from a test involving a set of 1000 simulations of peak detection. In other 
words, we fit a chunk of PSD by considering on one side a model of the peak with the background and 
on the other side a model of the background alone. The best model selected with the Bayesian model 
comparison encompassed in Eq. ^ is then telling us whether the peak is detected or not. 

The method presented in this work proves to be very efficient for this type of RG and will therefore 
be of great usefulness for the analysis of thousands of stars observed by Kepler . 
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Fig. 2. Linewidths of resolved dipole mixed modes for KIC 12008916 as a function of the frequency in the PSD. 



Fig. 3. Signature of the acoustic glitches for KIC 12008916 in the large frequency separation Av for the 
quadrupole and radial modes. 
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